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In the standard model (SM) the top-quark decay rate 
into a W boson and a down- type quark q (q = d, s,b) 
is proportional to |Vt g | 2 , the squared magnitude of the 
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element of the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix Under the assumption of a 3 x 3 unitary CKM 
matrix and using the existing constraints on Vt s and 
Vtd, the magnitude of the top-bottom quark coupling 
is |V t6 | = 0.99915+°;°™°5 and the top quark decays 

almost exclusively to Wb final states. Any significant de- 
viation from the expected value would imply new physics: 
an extra generation of quarks, non-SM top-quark produc- 
tion, or non-SM background to top-quark production. A 
direct measurement of the magnitude of the Vtb matrix 
element can be obtained from the single-top-quark pro- 
duction cross-section Q, which is proportional to |Vth| 2 . 
The value of \V t b\ can also be extracted from the decay 
rate of pair-produced top quarks. We define R as the 
ratio of the branching fractions 
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B(t -> Wb) 
B(t -> Wq) ' 



(1) 



Given the unitarity of the CKM matrix and assum- 
ing three generations, R is indirectly determined by the 
knowledge of V ts and Vtd, 
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and is derived to be 0.99830t^ooo9 0- A deviation from 
this prediction would be an indication of non-SM physics. 

This Letter reports the first CDF simultaneous mea- 
surement of R and top-quark-pair-production cross sec- 
tion o~ti performed on a data sample corresponding to an 
integrated luminosity of 8.7 fb _1 collected with the CDF 
II detector Q at the Fermilab Tevatron pp Collider at 
center of mass energy s/s = 1.96 TeV. The analysis uses 
events with a lepton and multiple jets in the final state, 
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where one W boson coming from tt production decays 
into a quark and an antiquark while the second W boson 
decays into a charged lepton (electron or muon) and a 
neutrino. 

CDF has performed several measurements of R dur- 
ing both Run I and Run II, combining the lepton+jets 
final state with the dilepton final state, where both W 
bosons decay into leptons. The most recent publication 
reported R = 1.12+°;^(stat) +°;^(syst), and R > 0.61 
at 95% confidence level (C.L.) using 162 pb _1 of inte- 
grated luminosity @. The DO Collaboration measured 
R = 0.90±0.04 (stat+syst) and R > 0.79 at 95% C.L. 
using data from 5.4 fb _1 of integrated luminosity, in the 
lepton+jets and dilepton final states combined. 

The CDF II detector [BJ consists of a charged-particle 
tracking system in a magnetic field of 1.4 T, segmented 
electromagnetic and hadronic calorimeters with a point- 
ing geometry, and muon detectors. A silicon microstrip 
detector provides determination of charged-particle tra- 
jectories (tracking) over the radial range 1.5 to 28 cm, 
and is essential for the detection of displaced decay (sec- 
ondary) vertices. A three-level, online event-selection 
system (trigger) Q is used to select events with an elec- 
tron (muon) candidate in the central detector region 
(pseudorapidity \rj\ < 1.1) with E T (p T ) > 18 GeV 
(18 GeV/c), which form the data set for this analysis. 

The measurement of R is based on the determination 
of the number of 6-quark jets in tt events reconstructed 
in the lepton+jets final state. The lepton+jets signature 
consists of a high-pr charged electron (e) or muon (/x), 
large missing transverse energy $ T Q due to the unde- 
tected neutrino from the leptonic W decay, and at least 
three hadronic jets. Events containing muons are classi- 
fied according to the coverage of the detectors used for 
their identification as central, when \rj\ < 0.6, and for- 
ward, when 0.6 < |?7| < 1.0. Identification of jets coming 
from 6-quark fragmentation (6-jet tagging) is performed 
by the SECVTX algorithm, which is based on the recon- 
struction of secondary vertices displaced from the pri- 
mary pp interaction vertex and selects a sample enriched 
with jets originating from b quarks [lOj . The lepton+jets 
selection requirements are described in Ref. |10j . Briefly, 
the analysis requires the presence of one isolated lepton 
(e or /i) with Et greater than 20 GeV, $ T of at least 
20 GeV, and a minimum of three jets, reconstructed 
using a cone algorithm [ll| with radius AR = 0.4 in 
r] — (j) space @, within |?7| < 2.0. The jet Et-, after 
correcting for the calorimeter response has to ex- 
ceed 30, 25, and 20 GeV for the most-energetic, second- 
most-energetic, and any additional jet in the event, re- 
spectively. The I^-boson transverse mass @ is required 
to be greater than 20 GeV/c 2 . Events with one or two 
identified 6-jets are selected (1 &-tag and 2 &-tag events, 
respectively). 

The background processes include VF-boson produc- 
tion in association with heavy-flavor jets (Wbb, Wcc, 
Wc), VF-boson production in association with light-flavor 
jets that are incorrectly identified as &-jets ("mistags"), 



quantum chromodynamics multijet ("QCD") events con- 
taining misreconstructed or real leptons or incorrectly- 
measured $ T , diboson events (WW, WZ, ZZ), single- 
top-quark production, and Z+jets events. 

We divide the selected sample into subsets according 
to the type of lepton, number of jets in the final state, 
and number of identified b jets (one or two). As explained 
in more detail below, we derive an expected event yield 
for each category. We then maximize the likelihood for 
observing the events found in each category by varying 
two fit parameters, R and the top-quark-pair-production 
cross section u t i- 

The events are modeled using the pythia 12] Monte 
Carlo (MC) generator with top-quark mass mt = 172.5 
GeV/c 2 . We estimate the backgrounds with a collec- 
tion of data-driven and simulation techniques described 
in detail in Ref. [l(|. The QCD background is modeled 
using data control samples |13j |. Mistags are estimated 
using a matrix (the mistag matrix) calculated in control 
samples and parametrized as a function of jet and event 
characteristics [Hj]. Diboson processes are simulated us- 
ing pythia; single-top-quark production is simulated by 



powheg [14| , while the parton shower and fragmentation 
is provided by pythia. The ALPGEN [ll| generator, with 
pythia supplying the parton shower and fragmentation, 
is used to model the VF+jets and Z+jets backgrounds. 
A GEANT-based simulation is used to model the response 
of the CDF II detector |16| . The cross sections used for 
background normalization can be found in Ref. [lH ] . Ta- 
ble U shows the expected sample composition for all final 
states, after summing over lepton categories, assuming 
R = 1 and cxti = 7.04 ± 0.49 pb [l7|- 

The number of fe-tagged events is the most sensitive 
quantity to possible values of R different from one: the 
smaller R, the smaller the probability to have a &-jet in a 
top-quark-pair event. Hence, the fraction of events with 
one or two tags is expected to decrease with decreasing 
R. In general, the tt production cross-section measured 
by CDF in the lepton+jets sample assumes R=l. In 
order to avoid any bias due to this premise, we measure 
simultaneously R and the production cross-section, since 
the measurement of the latter is affected by the sum of 
events in the different tag bins. 

To perform the fit, we first divide the sample into 18 
independent subsamples, organized by type of lepton, 
number of jets in the event (3, 4, > 5), and number of 
identified 6-jcts. The expected number of events in each 
subsample, /x^ , is given by the following expression: 



i. j 



^ e evt a ttel ag (R) 



(3) 



where [i^ is the expected number of tt events and Ng 3 
the expected number of background events. The i and 
j indices indicate the ith jet bin with one or two identi- 
fied 6-jets (jet-tag category) and the jth lepton category, 
respectively. Jzf J is the integrated luminosity, e^ 3 t in- 
cludes the trigger and lepton identification efficiencies, 
and e\ ag (R) is the event-tagging efficiency, i.e., the ef- 
ficiency for tagging at least one jet in an event. In an 
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TABLE I: Number of expected and observed events in lepton+jets data corresponding to 8.7 fb 1 of integrated luminosity. 
Uncertainties include statistical and systematic contributions. 







1 6-ta; 


:>' 




2 6-tags 




Process 


3 Jets 


4 Jets 


>5 Jets 


3 Jets 


4 Jets 


>5 Jets 


tt 


800 ± 67 


777 ± 64 


260 ± 21 


216 ± 30 


271 ± 36 


97 ± 13 


W+bb 


291 ± 118 


74 ± 30 


17 ± 7 


48 ± 20 


14 ± 6 


4 ± 2 


W+cc 


167 ± 68 


47 ± 20 


12 ± 5 


5 ± 2 


2 ± 1 


0.8 ± 0.4 


W+c 


87 ± 35 


17 ± 7 


4 ± 2 


3 ± 1 


0.8 ± 0.4 


0.2 ± 0.1 


Single top 


78 ± 7 


17 ± 2 


3.6 ± 0.3 


18 ± 3 


4.7 ± 0.7 


1.1 ± 0.2 


Diboson 


45 ± 5 


11 ± 1 


3.1 ± 0.3 


3.1 ± 0.5 


0.9 ± 0.1 


0.30 ± 0.05 


Z+jets 


32 ± 3 


9.1 ± 0.9 


2.4 ± 0.2 


2.1 ± 0.2 


0.77 ± 0.08 


0.29 ± 0.03 


Mistags 


303 ± 42 


74 ± 14 


17 ± 6 


5 ± 1 


1.7 ± 0.4 


0.6 ± 0.2 


QCD 


125 ± 50 


35 ± 29 


10 ± 9 


6 ± 3 


0.1 ± 1.5 


0.1 ± 1.5 


Total prediction 


1928 ± 243 


1061 ± 93 


330 ± 28 


306 ± 40 


296 ± 38 


104 ± 13 


Observed 


1844 


1088 


339 


275 


273 


126 



ideal case without background and assuming a 6-tagging 
efficiency equal one for jets originating from b quarks and 
zero for jets originating from non-6 quarks, the number 
of expected events with two tags is proportional to R 2 , 
while the number of expected events with one tag is pro- 
portional to 2R(\ — R). The estimates for the background 
processes are calculated with various values of R. The 
differences with respect to the estimates obtained with 
R = 1 arc found to be negligible. 

The event-tagging efficiencies are calculated in tt MC 
samples, using the probability to tag a jet as a &-jet ac- 
cording to the SECVTX algorithm, on a jet-by-jet basis. 
For jets originated from b and c quarks, the &-jet tagging 
efficiencies are corrected for differences between data and 
MC using a scale factor SF = 0.96 ± 0.05 [H|. For jets 
originated from light-flavor quarks, the probability to tag 
them as 6-jets is obtained using the mistag matrix. 

In general, t\ ag is calculated from the event probability 
to tag the mth event in MC processes with possible b- 
quark final states. For an event with n generic jets, the 
probability to have one (Eq. (@|) or two (Eq. ([5])) tagged 
&-jets becomes 

n I n \ 

p i n -t ag =Y.pT\ n i 1 -^ 9 )) (4) 

q— 1 \ r— l.r^q J 

7i— 1 / n I n \ \ 

c=e< 9 E^ as n (i-rf*) > 

q=l \r>q \s=l,s=£q,s=£r J J 

(5) 

where p*" 9 is either the probability to tag the gth jet, mul- 
tiplied by the SF, for jets where the heavy-flavor quark 
is found inside the jet cone; or the mistag probability for 
jets matched to a light-flavor hadron, calculated using 
the mistag matrix. 

Finally, we use the P™ tag as an event weight to calcu- 
late the event-tagging efficiency t\ ag for each subsample 
with I tags and n jets by summing the P™ tag weights 
over all of the pretagged events. 



The MC sample employed for the ti signal modeling 
is generated assuming |Vtf>| = 1 so it cannot be used di- 
rectly to calculate e\ ag as a function of R through the 
algorithm described above. Instead, in the MC sample 
we assign a random number p, in the interval [0, 1] to 
every jet that is matched at the parton level to a b quark 
from i-quark decay. If p, < R we consider this jet as gen- 
uinely originated by a b quark and use the tag probability 
multiplied by the SF as in Eq. (H|) and otherwise this 
jet is regarded as a light-flavor jet. This simulates a con- 
figuration in which a b quark produced in the top decay 
is a real b only R fraction of the time while (1 — R) frac- 
tion of the time it is treated as a light-flavor quark and it 
is weighted by the mistag probability. This probabilistic 
approach allows the calculation of background and sig- 
nal sample composition for any value of R. This method 
reproduces exactly the standard calculation in the case 
of R = 1, simulates t — s- Wq for R = 0, and allows a 
calculation of e l ta (R) in each tag subsample and in each 
jet bin. Figure [Tf shows the comparison of observed and 
expected events assuming R equal to 1, 0.5, and 0.1 in 
the various jet multiplicities and number of ^-tagged jets 
in the final state. 

In order to compare the prediction to the observed 
data in the 18 subsamples we use a likelihood function. 
We fit the observed event yields in each class of events 
to determine simultaneously R and a t i- The likelihood 
function is: 

c = n & {rfdpiR, °t-u x a )\Ki) n g { Xa % i) , (6) 

i,j a- 

where ^ (^fi l ^ p (R, a t i, x a )\N^ s ^j is the Poisson probabil- 
ity to observe N^f events assuming the expected mean 
Meip (given by Eq. ([3])), the index i indicates the jet-tag 
category, and the index j runs over the different lepton 
categories. The estimates of the nuisance parameters 
constrained to their central values and normal- 
ized to their uncertainties using Gaussian distributions 
G (x a |0, 1) centered at zero with unit variance. This pro- 
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FIG. 1: Observed events for the analysis final states after sum- 
ming over lepton categories, compared to expected events for 
different values of R. For the tt normalization the theoretical 
value a t i = 7.04 ± 0.49 pb is used. 



ccdurc takes into account correlations among channels by 
using same parameters for common sources of systematic 
uncertainties and allowing variations of each parameter 
with respect to its central value. 

We perform the minimization of the negative loga- 
rithm of the likelihood —2 log (£), using the MINUIT pack- 
age [jjj]. We analitically extend e\ ag (R) beyond R = 1 
during the fitting procedure, constraining each individ- 
ual e\ ag {R) to be greater than zero and their sum to be 
< 1. We simultaneously fit R and <7 t j, which are the free 
parameters of the likelihood. We update the calculation 
of background yields using the value of a t i determined 
by the fit and iterate the previous steps until the pro- 
cedure converges. No dependence on the starting point 
was observed in the results of the iterative procedure. 

The uncertainty determined by the fit comprises the 
statistical contribution; the systematic contribution on 
event-tagging efficiency, due to the systematic uncer- 
tainty on SF and the mistag matrix; the event selection 
efficiency, due to the lepton-identification scale-factor 
and the trigger efficiency; the background normaliza- 
tions, including the heavy-flavor fractions; corrections for 
differences between MC and data heavy-flavor yields; and 
the luminosity [l3j]. We include separately the contribu- 
tions due to the uncertainty on the jet-energy scale, ef- 
fect of initial- and final-state radiation in the simulation 
(ISR/FSR), event-generator dependences, and top-quark 
mass. The impact of the jet-energy scale uncertainty is 
estimated by varying the energy of all jets in the MC 
samples by ±1 standard deviation with respect to the 
central value for both signal and backgrounds and by re- 
peating the iterative fits. The uncertainty arising from 
the choice of the MC generator is evaluated by repeating 
the analysis using a tt sample generated by herwig [2(| ■ 
The ISR/FSR uncertainty is evaluated by using tt MC 
samples generated with enhanced or suppressed radia- 
tion relative to the default configuration. The theoret- 



TABLE II: Systematic uncertainties on the measurement of R 
and a t i obtained from simulated experiments. "Others" indi- 
cates the squared sum of minor systematic uncertainties. All 
systematic uncertanties are assumed to be fully uncorrelated. 
The statistical uncertainty is shown as well. 



Source 


+SR 


— on 


J-A/T - 

-r 0(J tt 


oa tt 








(Pb) 


(Pb) 


^-tagging 


0.078 


-0.073 


0.06 


-0.03 


Background normalization 0.056 


-0.052 


0.78 


-0.66 


Jet-energy scale 


0.016 


-0.019 


0.46 


-0.41 


ISR/FSR 


0.006 


-0.006 


0.22 


-0.21 


Luminosity 


0.001 


-0.002 


0.44 


-0.39 


Top-quark mass 


0.001 


-0.000 


0.33 


-0.32 


Others 


0.005 


-0.006 


0.17 


-0.15 


Total syst. uncert. 


0.088 


-0.081 


1.04 


-0.92 


Statistical 


0.043 


-0.043 


0.29 


-0.29 



ical value of the top-quark-production cross section de- 
pends on top-quark mass [2l| . The recursive fit of tr t j is 
expected to reduce the impact of this systematic uncer- 
tainty. In order to check this assumption, we repeat the 
measurement using two different MC samples for the tt 
signal, simulated with nit = 170 GeV/c 2 and m t = 175 
GeV/c 2 , respectively. Central values and uncertainties 
on those systematic effects are included in the likelihood 
as nuisance parameters. 

As a consistency check, the effect of each source of 
systematic uncertainty is estimated via simulated ex- 
periments. For each source we generate a set of simu- 
lated experiments with the same prescription but with 
the nuisance parameter x a , relative to the systematic ef- 
fect under study, shifted by one standard deviation from 
its nominal value. We determine the effect of changing 
each source of systematic uncertainty as the change in 
the mean of the distributions of R and cr t j. Table [TT] lists 
the various systematic uncertainties assumed as fully un- 
correlated. 

The final results are summarized in Table IIII1 Figure [2] 
shows the two-dimensional likelihood contour in the (i?, 
a t t) plane, for the fit including statistical and system- 
atic uncertainties. The best-fit value is indicated with 
a "X" and can be compared to the theoretical SM pre- 
diction at the next-to-leading (NLO) order expansion in 
the strong- interaction coupling constant [l7[ • The results 
are in agreement with the theoretical prediction to within 
one standard deviation. 

To determine the credibility level limit on R we fol- 
low a Baycsian statistical approach. Since R is bounded 
to be in the interval [0,1], the prior probability density 
is chosen to be zero outside these R boundaries, while 
we consider all physical values equally probable. To ob- 
tain the posterior distribution for R, we integrate over all 
nuisance parameters using non-negative normal distribu- 
tions as priors. We also integrate over a t t with the only 
constraint to be positive defined. The Baycsian lower lim- 
its at 68% and 95% credibility levels are shown in Fig. [3] 
and yield R > 0.785 at 95% c.l. From Eq. (TT]) we ex- 
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FIG. 2: The fit results for the simultaneous measurement 
of R and er t j. The X-cross corresponds to the maximum of 
the likelihood; the point with error bar to the NLO cross 
section calculation. The two-dimensional confidence regions 
are shown as well. 



tract a measurement of V tb - Assuming three generations 
of quarks and given the unitarity of the CKM matrix, we 
have \V td \ 2 + \V ts \ 2 + \V tb \ 2 = 1, leading to R = \V tb \ 2 . 
From the fit results we obtain |Vtb| = 0.97 ± 0.05 and 
\V tb \ > 0.89 at 95% credibility levels. 



TABLE III: Measured values resulting from the simultaneous 
likelihood fit of a t j and R. The uncertainties on 7? and o t t 
correspond to a variation of one unit of —2 log (£). The corre- 
lation parameter is p — —0.434. The magnitude of the CKM 
matrix element |Vtb| is derived from R = |Vti,| 2 . Lower limits 
at different credibility levels (c.l.) are obtained by integration 
of the posterior probability distribution. 



Parameter 


Value 


Lower limit 


Lower limit 




(stat+syst) 


68% c.l. 


95% c.l. 


0tt (pb) 


7.5 ± 1.0 






R 


0.94 ±0.09 


0.876 


0.785 


\Vtb\ 


0.97 ±0.05 


0.936 


0.886 



In summary, we present the simultaneous measurement 
of R = 0.94±0.09 and o t t = 7.5±1.0 pb with a correlation 
p = -0.434, and the determination of \V t b\ = 0.97±0.05. 
The results for R and |Vt&| are the most precise deter- 
mination obtained by CDF and are in agreement with 



the standard model [|J , with the previous CDF measure- 
ments ly], with the latest measurement of R performed 
by DO [7j , and with the direct measurement of single-top- 
quark production cross section performed by LHC [22| 
and Tevatron [H, [23[ experiments. 
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